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LOWER CRITICAL FIELD MEASUREMENTS IN 
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The temperature dependence of the lower critical field in 
YBa 2 Cu 3 0 6+x single crystals has been determined by magnetization 
measurements with the applied field parallel and perpendicular to the 
c-axis. Results were compared with data from the literature and 
fitted to Ginzberg-Landau equations by assuming a linear dependence 
of the parameter /c on temperature. A value of 7 i 2 kOe was 
estimated for the thermodynamic critical field at T = 0 by comparison 
of calculated H c2 values with experimental data from the literature. 


INTRODUCTION 

The lower critical field (H cl ) of the 
high temperature superconductor 
YBa 2 Cu 3 0 6 + x (YBCO) is an intrinsic 
material property which depends upon 
temperature and crystallographic 
orientation. Accurate measurement of H c 1 
is complicated by flux pinning and edge 
effects (as illustrated in Fig. 1) and by 
uncertainty in the demagnetizing factor. 
Consequently, early reported H c x values 
in single crystals 1 " 5 were up to an 
order of magnitude larger than later 
values . 6 ' 1 2 

In this study, H c 1 measurements were made 
on twinned and detwinned crystals and 
the data were compared with results of 
previous investigations for orientations 
with the applied field parallel and 
perpendicular to the c-axis of the 
crystal. The experimental data were 
fitted to Ginzberg-Landau equations for 
the dependence of H c x on temperature and 
upper critical field (H c2 ) values were 
calculated for the two orientations of 
interest . 

EXPERIMENTAL PROCEDURE 

The YBCO crystals used in the present 
study (Fig. 2) were grown from Y-Ba-Cu-0 


melts 13 and subsequently annealed in 
oxygen gas at 420 °C for 80 h to obtain 
superconducting transition temperatures 
T c > 90 K. (Thus, the oxygen content 
6+x > 6.85 14 ). Two crystals were 
selected for measurement. The first 
crystal (AN3-5) exhibited characteristic 
(110) twin planes and was nearly cubic 
with dimensions 120 x 135 x 120 /im 3 (c- 
dimension = 120 /^m) . Due to the cubic 
morphology, the demagnetizing factors 
were nearly identical in all three 
dimensions. The second crystal (AN9-5) 
was fully detwinned via a thermo - 
mechanical process developed in our 
laboratory 15 and had dimensions a x b x 
c = 200 x 250 x 100 fim 3 . 

Magnetic measurements were made using a 
superconducting quantum interference 
device (SQUID) magnetometer, with the c- 
axis of the crystal aligned either 
perpendicular or parallel to the applied 
field H. The crystal was first cooled 
in zero field to a predetermined 
temperature and the magnetization was 
then measured as the applied field was 
increased to a value in excess of H c 1 . 
For temperatures greater than about 60 
K, sharp breaks from linearity in the M 
vs . H curves were observed for both 
crystals, making the estimation of H c x 
relatively precise. Below 60 K, the 
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1 ‘ Hypothetical magnetization curves after cooling in zero field for a sample: 
(a) at equilibrium (i.e., an ellipsoid with no flux pinning); (b) with some pinning 
and edge effects; (c) with strong pinning; and (d) with perfect diamagnetic character 
(i.e. , magnetization is proportional to the applied field after correction for the 
demagnetizing factor, D) . For the equilibrium case (a), there is a well-defined, 
sharp break at H c 1 , which allows for an accurate determination of H c x . Pinning and 
edge effects (b and c) make it difficult to estimate the true H c 1 . In case (b) , the 
observed onset occurs at applied fields below H c 1 due to flux penetration at the 
sharp edges and corners of the crystal. In case (c) , pinning causes a gradual 
departure from linearity, making the estimate of H c ^ less accurate. 



Fig. 2. Scanning electron 
micrograph of as -grown YBCO 
single crystals. The smallest 
dimension of a crystal generally 
lies along the c-axis of the 
unit cell. 


250 


ORIGINAL PAGE 

BLACK AND WHITE PHOTOGRAPH 



departure from linearity was more 
gradual and H c x was estimated from the 
initial point of departure from 
linearity. 

RESULTS AND DISCUSSION 

Temperature -dependent H c x data from the 
present study and previous investiga- 
tions 4 ' 6 " 8 r 1 0 " 1 2 for H || c and H X c are 
presented in Fig. 3. The earliest 
reported H c x values 1 * 3 ' 5 were 
erroneously high due to difficulties in 
defining H c 1 and are not included in the 
two plots. The curves shown in each plot 
were obtained by fitting all displayed 
data points using the Ginzberg-Landau 
equation 1 6 

H cl = H c (In /c + 0.08)/72 k, (1) 

where H c is the thermodynamic critical 
field at temperature T as given by 

H c = H c0 (1 - t 2 ). (2) 

Here H c0 is H c at T = 0, t is the reduced 
temperature T/T c and k is the Ginzberg- 
Landau parameter (the ratio of the 
penetration depth A to the coherence 
length |). Data for both H || c and H l c 
can be well-fitted by assuming that k 
varies linearly with temperature: 

k = a + bt . (3) 

Good fits of Eq. (1) may be obtained for 
a wide range of H c 0 values, leading to a 
wide range of values for the parameters a 
and b in Eq. (3). The range of 
permissible H c 0 values is limited 
considerably by requiring that H c2 values 
calculated from a second Ginzberg-Landau 
equation 

H c 2 = J2 K ti c (4) 

be in reasonable agreement with 
experimental H c 2 data from the litera- 
ture. This comparative analysis yields a 
value of 7 kOe for H c 0 , k = 100 + 85t for 
H l c and k = 22 + 22t for H || c. The 


calculated curves and data for H c2 are 
compared in Fig. 4. Considering the 
large uncertainties in the experimental 
H c2 data and the obvious differences in 
T c for samples from different studies, 
agreement of the calculated curves with 
the data is a matter of judgement. The 
maxima seen in the calculated H c2 curves 
are unphysical, indicating that the 
linear form used for k should be 
modified, e.g., by the addition of a 
quadratic term in t. However, the large 
uncertainty in the H c0 value used here 
(7 ± 2 kOe) does not justify such an 
additional term. Considering the large 
variations in the experimental H c 2 data, 
our H c0 value is in reasonable agreement 
with the value of 10 kOe estimated by 
Worthington et al. 6 The resulting 
uncertainties in k values calculated 
from our equations are also of order 
±30%. 

Our H cl data for the detwinned crystal 
AN9-5 shown in Fig. 3a (H || c) are in 
good agreement with the data of Krusin- 
Elbaum et al. 10 for twinned crystals. 
This result indicates that twin bounda- 
ries have only a small effect on H c ^ , as 
noted in our earlier investigation. 17 
Our H cl data for the twinned crystal 
AN3-5 shown in Fig. 3b (H 1 c) are in 
reasonable agreement with the data from 
previous investigations / . 6-8, 10-12 
Anisotropy in H cl (H cl || c/H c x 1 c) as 
calculated from the two curves in Fig. 3 
was 3.1 ± 0.1 for 10 K < T < 80 K. 

The calculated H c2 values for H || c 
(Fig. 4a) are in reasonable agreement 
with the data of Welp et al. 18 near the 
superconducting transition and follow 
the general trend of the data of lye et 
al. 19 and Worthington et al. 2 . For Hlc 
(Fig. 4b) the calculated values at high 
temperature show good agreement with the 
experimental data of Gallagher et al. 3 
and Welp et al. 18 Anisotropy in H c2 
(H c 2 1 c/H c 2 || c) as calculated from the 
two curves in Fig. 4 was 4.3 ± 0.2 for 
10 K < T < 80 K. 
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- — l-i- Temperature -dependent H c 1 data for (a) H || c (strong pinning) and (b) H ± c 
(weak pinning). Our data are for (a) the detwinned crystal AN9-5 and (b) the twinned 
crystal AN3-5. The curves were generated by fitting the data to the Ginzberg- Landau 
equation for H c x and temperature -dependent k equations given by (a) k - 22 + 22t and 
(b) k = 100 + 85t (t = T/T ). ' ' 
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Fig. 4. 


Temperature - dependent H t2 curve (solid line) calculated from the^lnzberg- 


Landau equation for H o2 and the temperature -dependent * equation^forOi) H II c and 


1UI1 IUL n_ 9 CL I ILL die e— IT - - . rpA 

<b) H 1 c. Experimental H t2 data from the literature are shorn for comparison The 
extrapolations to lower temperatures (T < 60 K) are unreliable and the maxima 
curves probably do not exist. 




CONCLUSIONS 

Temperature -dependent H c 1 results were 
determined from magnetization measure- 
ments on de twinned and twinned single 
crystals of YBCO for H || c and H 1 c. 

The results from the present study and 
previous investigations for each 
orientation were fitted to Ginzberg- 
Landau equations assuming a linear 
temperature dependence for the parameter 

k. H c2 values calculated from the 
Ginzberg -Landau equation and the 
temperature -dependent k relations were in 
reasonable agreement with experimental 
H c2 data from the literature near the 
superconducting transition temperature. 
Values of H c0 = 7 ± 2 kOe , k = 100 + 85t 
for H l c, and k = 22 4- 22t for H || c 
were estimated from the analysis. 
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